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For lattice models with compact field integration (nonlinear sigma models over compact manifolds 
and gauge theories with compact groups) and satisfying some discrete symmetry, the change of sign 
of the bare coupling g$ at zero results in a mere discontinuity in the average energy rather than 
the catastrophic instability occurring in theories with integration over arbitrarily large fields. This 
indicates that the large order of perturbative series and the non-perturbative contributions should 
have unexpected features. Using the large- N limit of 2-dimensional nonlinear O(N) sigma model, 
we discuss the complex singularities of the average energy for complex 't Hooft coupling A* = <?oJV. 
A striking difference with the usual situation is the absence of cut along the negative real axis. We 
show that the zeros of the partition function can only be inside a clover shape region of the complex 
A' plane. We calculate the density of states and use the result to verify numerically the statement 
about the zeros. We propose dispersive representations of the derivatives of the average energy for 
an approximate expression of the discontinuity. The discontinuity is purely non-perturbative and 
contributions at small negative coupling in one dispersive representation are essential to guarantee 
that the derivatives become exponentially small when A* — > + . We discuss the implications for 
gauge theories. 



PACS numbers: 11.15.-q, 11. 15. Ha, 11.15.Me, 12.38.Cy 
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I. INTRODUCTION 

The lattice formulation of quantum chromodynam- 
ics provides a non-perturbative ultraviolet regularization 
and is a widely accepted model for strong interactions. 
For pure gauge SU(2) and SU(3) with a Wilson action, 
the absence of phase transition for any real and positive 
values of g% suggests that it is possible to use a weak cou- 
pling expansion in the bare coupling which is valid a 
short distance, to describe the large distance behavior of 
the theory. However, possible limitations of the valid- 
ity of the weak coupling expansion have been raised by 
Dyson He argued that if in quantum electrodynam- 
ics e 2 is changed into — e 2 , the vacuum becomes unstable 
and that consequently the radius of convergence of the 
expansion in powers of e 2 should be zero. In that article, 
Dyson says that "the argument [...] is lacking in math- 
ematical rigor [...] it is intended [...] to serve as a basis 
for further discussions". The idea was followed benefi- 
cially by Bender and Wu 0, [1[ and many others [H to 
show that the large order behavior of perturbative series 
can be estimated by semi-classical calculations at small 
negative coupling. 

The factorial growth of perturbative coefficients is due 
to the large field contributions to the path-integral [!, @ . 
In lattice gauge theory with compact groups, there is a 
build-in large field cutoff and the theory is well defined 
at negative g$. For pure gauge theory with a SU{2N) 
gauge group and lattices with even number of sites in 
each direction, there is an exact discrete symmetry [7| 



* yannick-meuriceOuiowa. edu 



that relates the partition function at opposite values of 
cjq. This symmetry implies that the average plaquette 
jumps suddenly from to 2 as g 2 goes from very small 
positive values to very small negative values. So, as <? 2 
changes sign, there is a change in vacuum rather than a 
loss of vacuum. 

The idea that Dyson argument needs to be revisited for 
lattice models with compact gauge groups is also sup- 
ported by the analysis of the expansion of the average 
plaquette in power of On. Existing series for SU(3) up 
to order 10 [8| and 16 Q suggest a power growth rather 
than a factorial growth. These perturbative series are 
constructed by adding the tails of integration that are 
absent because of the compactness of the group. In the 
case of a single plaquette [l(|, adding the tail of inte- 
gration leads to a factor ial g rowth. At finite volume, the 
studies of Refs. [ll[ and fl2j] suggest that factorial behav- 
ior may show up at an order proportional to the volume. 
The large order behavior of the perturbative series is re- 
lated to the zeros of the partition function in the complex 
l/g 2 plane [H, [3 (Fisher's zeros). This is clearly a dif- 
ficult problem and it would be useful to understand the 
connection for lattice models that have the same features 
but where calculations are easier. 

In this article, we discuss Dyson instability for nonlin- 
ear O(N) lattice sigma models. These models are defined 
in Sec. |TT] where we show that their partition functions 
have the same property as SU(2N) gauge theories under 
the transformation g^ —> —g$. In the large- N limit, we 
can use a saddle point approximation. This provides an 
equation relating the mass gap to the 't Hooft coupling 
A f = g^N. The complex singularities of this map and 
its inverse are discussed in Sec. Mil and ITVl respectively. 
An important feature is the absence of cut along the real 
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negative axis in the A* plane. 

The average energy and the density of states are cal- 
culated in Sec. [V] for large- TV and compared with weak 
and strong coupling expansions. In Sec. IVI1 we discuss 
Fisher's zeros. We show that they are related to the poles 
of the average energy and that they can only appear in 
a clover shape region of the complex A* plane. We use 
the density of states to numerically verify this statement. 
In Sec. IVII1 dispersive methods are proposed to repre- 
sent the derivatives of the average energy in the limit 
where A* — ► + . These representations are characterized 
by large contributions canceling each other. 

The article mentions other questions that would be 
worth addressing in more detail. One is the meaning 
of the finite radius of convergence for the linear sigma 
models. This seems to contradict Dyson's argument and 
may be relevant to understand the question of conformal 
fixed points suggested by Polyakov 15]. The other is 
the volume dependence of the non-perturbative part of 
the average energy and of the locations of Fisher's zeros. 
These questions are important for numerical calculations 
in lattice gauge theory. 



II. THE MODEL 



A. Basic definitions 

In this article, we consider the O(N) nonlinear sigma 
model on a square lattice. Most of the results presented 
in this section can be formulated for arbitrary dimension 
D not specified until the next section. We call V = L D 
the number of sites. The lattice sites are denoted x and 
the scalar fields </> x are A-dimensional unit vectors. The 
partition function reads: 

z = cJ]J d N ^8{$a* - i)e-^o)mm] , (i) 



B. Negative coupling duality 

Unlike the linear sigma model, the partition function 
is well defined at negative coupling. The X are unit 
vectors, and consequently, —1 < </> x -0x+e < 1- This 
means that the energy (or Euclidean action ) E per link 
is bounded from above and below. For a D-dimensional 
hypercubic lattice with an even number of sites in each 
direction and periodic boundary conditions, we have 



z[-gl] = ■ 



(•5) 



This can be seen by changing variable (f> — ► — on sublat- 
tices with lattice spacing twice larger and such that they 
share exactly one site with each link of the original lat- 
tice. A similar relation can be proven for SU (2N) pure 
gauge theories on even lattices [7|- 

The argument extends to compact manifolds (for sigma 
models) and to compact groups (for gauge theories) pro- 
vided that it is possible to transform the integration vari- 
able into minus itself without affecting the integration 
measure. It should also be noted that if ^ 0, Z[g 2 ] = 
implies Z[—g 2 ] = 0. 

The average energy per unit of volume V and its 
derivatives provide important information about possi- 
ble phase transitions. We denote it £ =< E > /V. With 
our notations, 



£ = -(l/(VN))d\nZ/db 
The symmetry ([5]) implies the sum rule 

£{-gl)+£{gl) = 2D . 



(6) 



(7) 



Knowing that when g% — ► + , £ (g%) — > 0, the sum rule 
implies that if g 2 — > 0~, then £ (g 2 ,) — > 2D. In other 
words, there is a discontinuity in the average energy when 
g 2 changes sign. 



with 



E[{4>}} =-^(0x.0x+e-l) , 



(2) 



with e running over the D positively oriented unit lattice 
vectors. We introduce the 't Hooft coupling: 



A* = g 2 N , 



(3) 



that is kept constant as N becomes large. Its inverse is 
denoted 



b= 1/A* 



(4) 



The volume integration at each site is finite and 
equal to the hypersurface of a N — 1 dimensional 
sphere 2tt n/2 /T(N/2). With the normalization C = 
(T(N/2)/2tt n / 2 ) v , the partition function becomes 1 in 
the limit 6 = 0. This is the analog of having the Haar 
measure normalized to 1 in lattice gauge theory. 



C. The gap equation 

In the large- N limit, it possible to calculate the 
partition function in the saddle point approximation 
[l6l [l7l [T§ | . In the case of the nonlinear sigma model, 
one enforces the condition <fi x .(j) x = 1 using a Lagrange 
multiplier. The integration over 4> can then be done ex- 
actly. Varying the zero mode of the Lagrange multiplier, 
we obtain: 



b = 05 (M 2 



with 



<B(A/ 2 ) = n r 

„• I J — 7T 



w dki 



L =V-* 2 ^ 2(Ejii(l -cos^ + M 2 



(8) 



(9) 



M 2 is the saddle point value of the suitably rescaled La- 
grange multiplier and can be interpreted as the mass gap 
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or as the renormalized mass in cutoff units. At finite 
volume, the integral is replaced by a sum over momenta 
rii2ir/L. 

The saddle point equation is invariant under the simul- 
taneous changes: 



A* 
M 2 



-A* 
-M 2 



(10) 



4L> 



This can be seen by changing variables kj — > kj + tt for 
all j. Note that this change of variable sends the zero- 
momentum mode into the fastest oscillating one (that 
changes sign at every lattice site). 

It is interesting to compare the gap equation (JSJ) with 
its counterpart for the linear sigma model. In the linear 
case, the Lagrange multiplier is used to replace (f> x .<p x by 
a composite field. After suitable rescalings [l7| , this com- 
posite field is denoted X and the gap equation becomes 



X 
M 2 



<B(M 2 ) 
2U'(X) 



(11) 



for a rescaled bare potential U(X). In the case of a </> 4 
theory, we can choose U(X) — (m|/2)J + XX 2 , where 
m 2 B is the bare mass and A is kept constant when N 
becomes large. X can then be eliminated yielding: 



(M 2 -m 2 B )/\ = *8{M 2 ) 



(12) 



Note that despite the fact that the partition function is 
not well defined at negative A, the D — 2 gap equation 
has 2 real solutions for M 2 when A c < A < 0. At A c 
the two solutions merge. They disappear in the complex 
plane when A < A c . 

The comparison between the linear and nonlinear 
sigma models shows that there are important differences 
between the two cases. In the linear case, the situation is 
very similar to what is observed in other scalar models [4] 
(dispersion relations with a cut extending to — oo), but 
the nonlinear case is quite different as we now proceed to 
explain. 



III. THE GAP EQUATION AT COMPLEX 
COUPLING FOR D = 2 

In this section, we study the map <B(M 2 ) when M 2 
is varied in a cut complex plane. The cut is between -8 
and on the real axis. D = 2 is assumed through this 
section. From the previous section, we know that 



Q3(-8-M 2 ) 
5B(M 2 *) 



-Q3(M 2 ) 
<8(M 2 )* 



(13) 



and we only need to study one quadrant, for instance 
Re M 2 > -4 and Ira M 2 > . 
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FIG. 1: Complex values of A* = 1/93(M 2 ) when M 2 varies 
over horizontal (with a spacing 1) and vertical lines (with a 
spacing 0.5) centered about the cut in the complex M 2 plane. 



A. Large M 2 behavior 

The basic observation is that for large \M 2 | , we have an 
approximate one-to-one mapping since A* ~ M 2 . This is 
fundamentally different from what happens for the linear 
sigma model where A ~ (1/4)(M 2 ) 2 and the inverse map- 
ping requires a cut that can be taken along the negative 
axis. 

The approximation can be improved and we can cal- 
culate the strong coupling expansion 

1/A* = 1 /M 2 - 4/M 2 + 20/M 6 . . . (14) 



and its inverse 



1/M 2 = l/A* + 4/(A*) 2 + 12/(A*) 3 



(15) 



Numerical studies of the large order expansion in pow- 
ers of 1/M 2 show a clear evidence for a singularity at 
1/M 2 = -1/8. 

As we reduce |M 2 |, the image of a cartesian grid gets 
distorted and gaps open near the origin. This illustrated 
in Fig. [T] The approach of the cut is intricate and in- 
volves logarithmic divergences in 55 (M 2 ) that we now 
proceed to study. 



B. Logarithmic divergences of 23 (M 2 ) 

For D = 2, logarithmic divergences appear in 5B(Af 2 ) 
from region of integration where 



2(2 — cos(fci) — cos(fc2)) 



j4. k-^ 



(16) 
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FIG. 2: Real (blue online) and Imaginary (red online) part 
of <8(M 2 ) when M 2 varies over a horizontal line 0.01 above 
the cut in the complex M 2 plane. The imaginary part has a 
negative spike at -4. 



This only occurs when fc, = or it. The four cases are A 
= (for (0,0)), 4 (for (0,tt) or (tt,0)) and 8 (for (n,Tr)). 
Logarithmic divergences appear when M 2 approaches 0, 
-4 and -8. In order to give a first idea, we have plotted 
in Fig. [2] the real and imaginary part of 23 (M ) when 
M 2 runs over a line slightly above the cut. This figure 
suggests that the real part has singularities at and -8 
and that the imaginary part has singularities at -4. 

The leading coefficients of the logarithmic singularities 
can be estimated by using the continuum approximation 
(UHJ). For instance, if M 2 — ► + on the real axis, we ob- 
tain the familiar relation Q3(M 2 ) ~ -(l/47r)ln(A/ 2 ). For 
M ~ —4, we have two contributions and by applying 
the proper Wick rotations, we obtain an imaginary ex- 
pression that is twice larger in absolute value than the 
real part near for M 2 ~ 0. 

Further insight into the singularities can be obtained 
by introducing the spectral decomposition : 



with 
G(u) 



<B(M 2 ) = / duG{u) 



cPk 

(27T) 



2vr 2 



1 



M 2 



(17) 



2 S(u - 2(2 - cos(fci) - cos(fe))) (18) 



dc 



[1- \u/2-2 + c\ 



^/^^ 2 v / l-( ^i /2-2 + c) 2 



The symmetry (|10[) implies that 

G(8 - u) = G(u) . 




FIG. 3: G(u) (blue online) compared to l/(47r) (red online) 



Numerical values for G(u) are provided in Fig. [3J The 
estimate of the singularity of the real part of 03 (A/ 2 ) near 
M 2 = implies 



G(0) = 1/4tt , 



(20) 



which is confirmed numerically in Fig. [3] From the 
residue theorem, we find that the imaginary part has a 
discontinuity along the cut. For M 2 G [—8,0], 



Im(23(M 2 ± ie)) = T^G{-M 2 ) 



(21) 



From the analysis of the logarithmic singularity of the 
imaginary part of 23 (A/ 2 ) near M 2 — —4, this implies 
that for u ~ 4, 



G( U )^~(l/(27r 2 ))ln(| W -4|). 



(22) 



This result was confirmed by an analysis of the numerical 
values of G(u) near u = 4. Finally, using the fact that 
when u reaches 4, two pairs of inverse square root sin- 
gularities coalesce into two poles, it is possible to justify 
that the real part jumps suddenly from -1/4 to 1/4 as 
M 2 increases and crosses -4. 

By construction, the integration path in Eq. (fT7)l does 
not wrap around the pole at u — —A/ 2 . As in the case of 
the logarithm, we can introduce 



03 fe (M 2 ) = Q3(Af 2 ) + ik2irG{-M 2 ) 



(23) 



(19) 



which corresponds to have the contour of integration 
wrapping k times around the pole. 
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FIG. 4: Complex values taken by 25 (M 2 ) when M 2 varies over 
the complex plane (here on horizontal lines in the M 2 plane 
with spacing 0.1 (black, blue online) and 0.5 (gray, orange 
online)). Asymptotic limits are ±0.25 in both direction. 



C. 03 image of the cut M plane 

We are now in position to determine the image 05 (M 2 ) 
of the cut M 2 plane. The region of large \M 2 \ is mapped 
into a neighborhood of the origin. As we approach the 
cut, the image of lines of constant imaginary M 2 are 
mapped into hat-shaped curves shown in Fig. dJ In the 
limit of zero imaginary part, the curves become approx- 
imate hyperbolas with asymptotes on the boundary of a 
cross of width 0.5 centered at the origin. The asymp- 
totes correspond to the logarithmic singularities and can 
be read from Fig. [2] When M 2 has zero real part and 
a small positive imaginary part, the imaginary part of 
03(M 2 ) reaches -0.25 while the real part becomes large 
and positive. The other asymptote of the approximate 
hyperbola is reached by approaching -4 from above with 
a positive imaginary part, the real part of 03 (M 2 ) then 
reaches 0.25 while the imaginary part becomes very neg- 
ative. All the other cases can be obtained from the sym- 
metry (TT3|) . 



D. 



A* image of the cut M 2 plane 



We can now draw more precisely the empty region in 
the middle of Fig. [1] . Before constructing it, we state 
the final result: the image of the cut complex M 2 plane 
under the 1/03 (M 2 ) map is the complex plane minus a 
clover shape centered at the origin and with lobes ap- 
proximately bisecting the real and imaginary axis. The 
clover shape shape is visible in Fig. [5J 




FIG. 5: Complex values taken by A' when M varies over 
lines above and below the cut with ImM 2 — ±0.01 (black, 
blue online); the circles are the inverses of the asymptotic 
lines in Fig. U 



The boundary of the clover shape is the 1/z map of 
the 4 limiting approximate hyperbolas discussed in the 
previous subsection. First, we construct the image of the 
asymptotes. Their images are circles of radius 2 centered 
at ±2 and ±2i. Near the origin in the A* plane, the 
circles provide a good approximation of the boundary. 
As we move away from the origin, numerical values are 
necessary. This is illustrated in Fig. [5l 



IV. SINGULARITIES OF M 2 (A*) 



In Sec. IIII1 we have constructed a one-to-one map 
between the cut M 2 plane and the A* plane with a clover 
removed. The clover can be filled by adding the images 
03fc indexed by an integer k as defined in Eq. (f2"5)) . In this 
Sec, we construct the inverse mapping and discuss its 
singularities. We start with a simplified example where 
everything can be done analytically and then discuss the 
original problem. 



A. A simplified example 

We will start by constructing the inverse mapping for 
a simplified mapping 03 S i mp j.(M 2 ) for which G(u) in the 
decomposition (fl~7|) is constant and equal to 1/4-7T. This 
implies 



KaimplXM 2 ) = (l/47r)ln(l ± 8/M 2 ) 



(24) 
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Poles of M 2 simpl . (X) 

i — 1 — t — 1 — i — 1 — i — ■ — r 

aes * 




I I i I i I i I i I I 

-4 -2 2 4 

Re X 

FIG. 6: A finite number of poles and boundaries of domains 
described in the text for X S im P i.(M 2 ) in the M 2 plane. 

This modification preserves the logarithmic singularities 
at and -8 and the symmetry (jTOJ) . This example is 
probably closer to what we expect to find for 4D lattice 
gauge theories. The image of the cut complex M 2 plane 
(fundamental domain) is the horizontal strip bounded 
by the horizontal lines with imaginary part ±0.25. The 
inverse map can be calculated explicitly: 

M 2 mp ,(&) = 8/(e 4 - b -l) , (25) 

and is invariant under b — > b+ik/2, for any integer k. The 
translated values obtained from the fundamental domain, 
could be obtained directly by having the path of integra- 
tion in Eq. (|17[) wrapping k times around ~M 2 . The 
poles of M 2 impl (b) are located at b = 0, ±i/2, ±i, .... 

In the A* plane, the fundamental domain is the com- 
plex plane with two circles removed . The boundary of 
the fc-translated domains are circles of radius 2/(2|fc| + 1) 
and centered at ±i2/(2\k\ + 1). The poles are located 
at —i2/k. A finite number of boundaries and poles are 
shown in Fig. [51 For large |fc|, the poles and boundaries 
accumulate at the origin. We believe that some qualita- 
tive features of this picture are representative of what is 
encountered in lattice gauge theory. 

B. Qualitative features of M 2 (A') 

We can now describe qualitatively what happens when 
we restore the original features of G(u). The fundamen- 
tal domain in the b complex plane becomes strongly dis- 
torted when Reb is small as already shown in Fig. 0] 
The translated domains that were parallel strips in the 
simplified example are obtained by adding ik2nG(—M 2 ) 
as specified by Eq. (|2"3")1 . The reader can visualize the 
effect by combining Figs. [Hand [3] 



We can now discuss qualitatively the deformation of 
Fig. [S]that the logarithmic divergence of G(u) near u = 4 
imposes. In the b plane the strips are pulled at infinity 
when Reb becomes small. In the A* plane, the circles 
arc pushed toward the origin along the imaginary axis 
forming concentric clover shape figures. Importantly, all 
the poles have moved to the origin. 

V. THE DENSITY OF STATES 
A. Average energy 

From the definition of the average energy Eq. ([S]), we 
obtain that in the saddle point approximation 

£ = (l/2)(\ t ~M 2 ) . (26) 

Note that < £ < 4 and the range is AMndependent. 

When A* approaches on the positive real axis, M 2 
approaches zero like 8e~ 47r /\ We call -(1/2)M 2 the 
non-perturbative part of £. The perturbative expansion 
terminates at first order. For large A*, the leading terms 
cancel and can use the strong coupling expansion (fT5"]) to 
obtain 

£ = 2-2/X t ... (27) 

B. Saddle point calculation of n(S) 

We can use a spectral decomposition of the partition 
function for all possible energies: 

Z = J dEn(E)e- bNE , (28) 

with n(E), the density of state which can be defined as 

n{E) = / II - - E) (29) 

n(E) is non-zero only if < E < DV in D dimensions 
and this implies that Z is an analytical function in the 
entire b plane. Using 

S(E[{<p}} -E)= dae a{E[W] - E) . (30) 

and varying with respect to M 2 (introduced as before) 
and a, we obtain 

a = »(l/a-2£) 
M 2 = l/a-2£ , (31) 

with a and M 2 understood as functions of £. In these 
equations £ is the independent variable and they are 
equivalent to Eqs. JSj) and (|26| provided that we identify 
a and b = 1/A*. From these results, we obtain 

n(E) oc e VN ^ E ' v ^ , (32) 
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with the entropy density 

f(£) = -(l/2)log(a) (33) 

2 2 

-a/2) n f ^m2(E( i - cos (^))+ m2 ) 

The duality ([5]) implies 

f(£) = f(4-£). (34) 

Numerical values of f(£) are shown in Fig. [71 They were 
calculated using M 2 as a parameter first used to fix a 
and £ from (011) and then f{£) from Eq. ([531). 

If we rewrite the partition function in terms of the 
density of states, we obtain the usual thermodynamic 
relation: 

f{£) = b . (35) 

One can check that the saddle point equations imply a = 
b as expected. 

The behavior of /(£) can be approximated near the 
origin. Using £ ~ (1/2) A* and Eq. ([33]) . we obtain 

f(£) ~ (l/2)ln(f ) . (36) 

The behavior of f(£) near its maximum at £ — 2 can be 
approximated using Eqs. ([27]) and ([35|) . The result is 

/(£)^(-l/4)(£-2) 2 (37) 

These two approximations fit the numerical values quite 
well in their region of validity as shown in Fig. [7J 

As long as £ takes values in the real interval [0, 4], Eq. 
(|35| defines a one-to-one mapping between this interval 
and the entire real axis. It is clear from Fig. [Jj that in 
this interval f"{£) < and consequently the derivative of 
the mapping is never zero over the interval. This allows 
us to integrate over the quadratic fluctuations about the 
minimum. If an analytical continuation of f{£ ) can be 
constructed, it seems clear that poles and zeros of f"(£) 
will play an important role in the determination of the 
zeros of the partition function which is the subject of the 
next section. 



VI. ZEROS OF THE PARTITION FUNCTION 

There exists a simple relation between the poles of the 
average energy and the zeros of the partition function. If 
bo is a zero of Z of order q, then (dZ/db)/Z ~ q/(b — bo) 
for b ~ &o- If we now integrate over a closed contour C, 

/ db{dZ/db)/Z = i2n V n q {C) , (38) 

where n q {C) is the number of zeros of order q inside C 
Using Eq. ([2"6")h we obtain that in the large N limit, 

(4m)- 1 I db(M 2 - l/b) =^2<n q {C)/(VN) . (39) 



In this expression, M 2 is understood as a function of b by 
inverting Eq. ([5]). The second term has a pole at b = 0, 
but it is compensated by a pole in M 2 . This is due to 
the fact the b ~ 1/M 2 for small \b\. We now consider 
possible poles of M 2 for other values of b. We change 
variable to write 

j> dbM 2 = I dM 2 (db/dM 2 )M 2 (40) 
Jc Jc 

where C is the contour corresponding to C in the M 2 
plane and b a short notation for 25 (M 2 ). At finite volume, 
it is possible write 25(M 2 ) as a ratio of two polynomials. 
The zeros of the denominator can only be in the cut (the 
real interval [—8,0]). This property persists for arbitrar- 
ily large volume. Consequently, if the contour C' in the 
M 2 plane does not cross the cut, then there are no ze- 
ros of the partition function inside the corresponding C 
in the b- plane. We conclude that in the large- N limit, 
there are no Fisher's zero in the 25(M 2 ) image of the cut 
M 2 plane. This image has been constructed in Sec. IIIII 
and limited by four approximate hyperbolas with asymp- 
totes along a cross shaped figure. In the A* plane, this 
region becomes the complement of a clover shape figure 
(see Figs, g] and 03). 

The argument has been checked by numerical calcula- 
tions using methods similar to those used in lattice gauge 
theories [H, [l3| . We used spline interpolations from 400 
numerical values of f(£) in Eq. ([33j and NV = 100. In 
order to remove fast oscillations, we have subtracted the 
average value of the energy at b = 0.35 from £ in the ex- 
ponential. This does not affect the complex zeros. The 



Entropy density 




E/V 

FIG. 7: f(£), numerical (circle), first order strong coupling 
(parabola, green online) and first order weak coupling (red 
online). 
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NV=10 




.45 



FIG. 8: Fisher's zeros for NV = 100. Zeros of ReZ (small 
dots, blue online), zeros of ImZ (larger dots, red online). The 
solid line (blue online) is the image of a horizontal line slightly 
below the cut in the M 2 plane. 



results are shown in Fig. [5] for a rectangle with 0.1 on 
each side of b = 0.35 and Imb < 0.4 (beyond that new 
methods need to be developed to calculate the rapidly 
oscillating integrals). The Fisher's complex zeros are the 
three isolated points at which the zeros of the imaginary 
part of Z meet those of the real part. The image of a line 
slightly below the cut (ImM 2 =-0.05) is also shown. As 
predicted there are no Fisher's zeros below the image of 
this line. 

The construction points out a relationship between 
Fisher's zeros and the poles of db/dM 2 . This derivative 
also appears in the /3-function 



/3(A*) ee AdX l /dA = 2(X t ) 2 M 2 db/dM 2 



(41) 



assuming M 2 



j /A 2 with the renormalized mass 



kept fixed and A the UV cutoff. In the approximation 
where b = — (47r) lnAf , we recover the well-known con- 
tinuum result /3(A 4 ) = — (A*) 2 /(27r). Following the finite 
volume reasoning used above, we see that the poles are 
exactly cancelled by zeros of (A*) 2 . The zeros of the beta 
function can be interpreted as the singular points of the 
A*(M 2 ) map which have nontrivial locations at finite vol- 
ume. 



VII. TOWARD A DISPERSIVE APPROACH 

In quantum mechanics or linear scalar models with a 
A</> 4 interactions, it is common 0, H, 0] to use dispersion 
relations to estimate the large order behavior of pertur- 
bative series. Typically, one considers a quantity, that 
we will denote F(X), that is analytical in the cut plane 



with a cut going from -co to 0. The coefficients Fk of 
the perturbative series in powers of A are then expressed 
as 



r° 

Ffc = (l/7r)/ dX ImF(X)/X k+1 

J — oo 



(42) 



For large k, the integral is dominated by small nega- 
tive values of A and one can use semi-classical meth- 
ods to estimate ImF(X) in this regime. Typically, 
ImF(X) ~ A~ b e a / A which leads to a factorial growth 
F k ~ (~a)~ k T(b + k). 

We are interested in finding a dispersive representation 
of £(X t ) as expressed in the leading order Eq. (|2"5j) . In 
this approximation, the series terminates and what we 
would expect to learn from dispersive methods is that 
M 2 , which appears in the second term off, is zero to all 
order in A, when A* — ► + . 

In the rest of this section we will study dispersive ex- 
pressions for M 2 impl defined in Eq. (|2"S")) . This simplified 
expression of the mass gap does has a different behav- 
ior when A approaches along the imaginary axis but it 
preserves the basic symmetry (110[) and it probably has 
more resemblance with the gauge models. This choice 
also has the advantage that all the calculations can be 
done in terms of elementary functions and that the cor- 
rectness and accuracy of dispersive expressions can be 
checked easily. We define 



: m (e) = (121T)- 1 



<j> dX(X - 



e)- m - X M* impl XX) , 



(43) 



where the contour of integration runs counterclockwise 
along a circle of center e and radius smaller than e. In 
the rest of this section, e is real and strictly positive. By 
construction, m\9\ m (e) is the m — th derivative of M 2 impl 
evaluated at e. If e is not exactly zero, we can obtain 
analytical expressions, such as 



JKi(e) = 87r/(sinh 2 (27r/e)e 2 ) . 



(44) 



The 9t m (e) calculated that way are even functions of e 
for m odd and larger than 1 and odd functions for m 
even and larger than 2. However, since M 2 impl has a 
discontinuity at 0, 8(e) and derivatives should also be 
present if we wish to obtain expressions valid over the 
whole real axis. 

For < e << Air, 9i m (e) can be approximated by a 
sum of terms of the form e~ 47T ^ e e~ q , with q < 2m. These 
individual terms peak at e = 4n/q where their value grow 
like q\. When m increases, there is an intermediate region 
where many large terms cancel and it is difficult to eval- 
uate the function numerically. For instance, for m = 12, 
some numerical noise becomes visible for 0.3 < e < 0.8 
when 16 digit arithmetic is used. 

For m > 2, the contribution of the circle at infinity 
vanishes, and we can deform the contour in Eq. (|43[) into 
a clockwise contour encircling the poles on the imaginary 
axis. Due to the absence of cut or poles on negative real 
axis, the contribution of two line integrals running on 
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opposite directions along the negative axis cancel. This 
is contrast to what happens for quantum mechanics and 
linear scalar models, where there is a discontinuity. Using 
the residue theorem, we obtain 

JMe) = -(8/t) ]T k m -\2i - ek)-" 1 - 1 . (45) 

This expression provides reasonably accurate expressions 
with truncated sums when e is not too small. When e 
becomes small, the sum can be divided into two parts, 
one with |fe| < 2/e and the other with the rest of the 
terms. Each sum is of order e _m and the two large sums 
must cancel to yield an exponentially small results. 

Another option is to start with the deformed contour 
encircling the poles on the imaginary axis keeping it 
slightly outside of the figure 8 visible on Fig. [6] Chang- 
ing variable to M 2 , this contour is then mapped into a 
contour encircling the cut from -8 to (but not not con- 
taining M 2 impl (e)). With this contour, we have 

m m (e) = {2m)- 1 j dMH'b™- 1 ^ - be)-" 1 - 1 , (46) 

with b a short notation for ^Bsimpi. (M 2 ) given in Eq. (pM)) 
and b' its derivative. The contour can be constructed as 
a "barbell" figure with three part that we call A, B and 
C: two lines running along the cut in opposite directions 
(part A) and two circles going around -8 (part B) and 
(part C) and stopping when they meet the horizontal 
lines. We have monitored the contributions of the three 
parts for m < 12 and values of e < 2 and compared them 
to reliable numerical values of the derivatives of M 2 impl . 
For small, but not too small values of e, A and B pro- 
vide same sign contributions that dominate the integral. 
As e decreases, the A and B contributions become large 
and of opposite sign while C becomes negligible. As e 
keeps decreasing, the absolute values of A and B keep 
increasing until the accuracy of the integrals deteriorate. 
The situation is illustrated for m — 5 in Table HI Again, 
we are facing the challenge of having two large canceling 
quantities. We have tried to improve this situation by 
modifying the radius of the circles and the distance be- 
tween the horizontal lines, but the general features seem 
quite persistent. 

The lesson that can be learned from the second rep- 
resentation is that the contributions from small negative 
coupling, which becomes the integral around a cut circle 
around M 2 = —8 (part B) after the change of variable, is 
essential to compensate the large contributions from the 
cut. 



VIII. CONCLUSIONS 

In summary, for 2D nonlinear O(N) sigma models, the 
exact discontinuity of the average energy appears to be a 
purely non-perturbative phenomenon. At leading order 



e 


T =9t 6 (e) 


A/T 


B/T 


C/T 


2.0 


-0.00116908 


-0.166729 


0.536511 


0.630219 


1.9 


-0.00202112 


0.257416 


0.37977 


0.362814 


1.8 


-0.00296429 


0.441387 


0.320456 


0.238158 


1.7 


-0.00389598 


0.527325 


0.305554 


0.167121 


1.6 


-0.00465005 


0.552561 


0.325353 


0.122087 


1.5 


-0.0049979 


0.517248 


0.390845 


0.0919072 


1.4 


-0.00467868 


0.379432 


0.548854 


0.0717137 


1.3 


-0.00347855 


-0.0515184 


0.99067 


0.0608484 


1.2 


-0.0013698 


-2.53523 


3.45779 


0.0774331 


1.1 


0.00131813 


6.10486 


-5.07936 


-0.0255078 


1.0 


0.00380812 


3.56912 


-2.56953 


0.000410349 


0.9 


0.00507521 


3.93213 


-2.93389 


0.00175748 


0.8 


0.00449056 


6.30264 


-5.30373 


0.0010864 


0.7 


0.00257563 


16.7544 


-15.7548 


0.000454331 


0.6 


0.000812222 


93.4553 


-92.4554 


0.000131085 


0.5 


0.0000987024 


1576.99 


-1576.15 


0.0000229591 



TABLE I: Values of the contributions A, B and C described 
in the text in units of the total value T = 915(e) • 



in the 1 /N expansion the series terminates and no patho- 
logical behavior at negative coupling can be inferred from 
it. It would interesting to see if this feature persists for 
sub-leading corrections. 

In the large- N and large volume limit, the Fisher's ze- 
ros can only appear inside a clover shape in the complex 
A* plane. It is plausible that as the volume increases, the 
zeros become dense at the boundary of the clover shape 
and at the boundary of an infinite set of concentric clover 
shapes. If this is correct, the zeros closest to the real axis 
in the b plane appear at infinite Reb and Imb — ±i/4. 
Numerical studies at finite N and V should clarify the 
picture and could be used as a guide for the search of 
Fisher's zeros in gauge theories. It would also be inter- 
esting to consider the case D = 3, where there is a rich 
phase diagram [It], EH, H3| and where we expect the zeros 
to pinch the real axis in the infinite volume limit. 

A simplified form of the discontinuity was approached 
with dispersive methods. An important feature observed 
was that large contributions cancelled. Again, finite vol- 
ume studies may clarify the mechanism. If we extend the 
reasoning used for one plaquette [Til ] . non perturbative 
effects become important at an order k ~ j3V. Near this 
order, the coefficients become different if we integrate a 
perturbative expansion of the density of state from to 
DV or from to oo. This understanding is crucial if we 
want to modify the weak coupling expansion in order to 
include non-perturbative effects. 

We have noticed that the zeros of the (3 function are 
related to the singular points of the mapping b(M 2 ). At 
finite volume, these singular points have a nonzero imag- 
inary part. This should be seen as an encouragement to 
study complexified renormalization group flows as also 
suggested by other work [IH, [2l[ ■ 
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